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ABSTRACT 



We present a series of four simulations of Cold Dark Matter (CDM) and Cold + Hot 
Dark Matter (CHDM) cosmologies, which we analyze together in this and subsequent pa- 
pers. These dissipationless simulations were done using the Particle Mesh method with a 
512 3 mesh, corresponding to a resolution of approximately 200 kpc for an assumed Hub- 
ble parameter of Ho = 50 km s 1 Mpc -1 , and with approximately 17 million cold and (for 
CHDM) an additional 34 million hot particles. In this paper we discuss the power spectrum 
and correlation functions in real and redshift space, with comparisons to the CfA2 and IRAS 
redshift data, the pairwise velocity of galaxies in real space, and the distribution of hot and 
cold particles in CHDM simulations. We confirm that CHDM with cold/hot /baryon density 
ratios fi c /fij, /fi& = 0.6/0.3/0.1 is a good fit to a wide variety of present-epoch data, much 
better than CDM. In particular, with reasonable assumptions about identification of galaxies 
and biasing, we find that the power spectrum from our CHDM simulations agrees rather well 
with both the CfA2 and IRAS power specta in both the nonlinear and linear regimes. 

New variants of the CHDM scenario (e.g., with 20% of the mass in hot particles or with 
two massive neutrinos) predict a significantly larger rate of formation of galaxies at high 
redshift, which may be needed to explain some observational data. At the same time, the 
difference between the variants is rather small at z = 0. The results presented in this paper 
are interesting for two purposes: (i) For a rough comparison with other classes of models 
(like CDM and ACDM) at z = — indeed, we have used the simulations described here 
as a test bed for developing a number of new statistics for quantifying large scale structure 
and comparing it to observations, (ii) As a reference point for comparison between differ- 
ent variants of the CHDM model. In addition, we explain here how we modify the usual 
Zel'dovich approximation used to set up the initial conditions for both cold and hot particles 
in our simulations, taking into account that the growth rates of both kinds of fluctuations 
are different from the usual CDM case. 

Subject headings: cosmology: theory — dark matter — large scale structure of the universe 
— galaxies: formation — galaxies: clustering 
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1. INTRODUCTION 



Cold Dark Matter (CDM) was perhaps the simplest cosmological theory proposed in the 
1980s that had a chance to be true. Standard CDM is based on the hypotheses that the 
primordial fluctuations are gaussian and adiabatic with a Zel'dovich spectrum, as predicted 
by simple inflation models, and that the dark matter is "cold," i.e., preserves fluctuations 
on all cosmologically relevant scales. Standard CDM is also usually assumed to have a total 
density which is critical (0 = 1). Aside from the Hubble parameter Ho = 100/i km s _1 Mpc -1 
and the fraction of critical density in ordinary matter (baryons), 0&, the only free parameter 
in CDM is then the overall normalization of the spectrum. Since for = 1 the age of the 
universe to = (2/3)H ~ 1 = 6.52/i -1 Gy is shorter than the age of the oldest globular clusters 
unless h < 0.5, it is usual to assume for = 1 that h = 0.5, and we make that assumption 
in this paper. 

In the earliest work on CDM (Peebles 1982; Blumenthal et al. 1984), the spectrum was 
normalized so that erg = 1, where 



the normalization is also sometimes specified by quoting the linear bias parameter b. The first 
CDM N-body simulations (Davis et al. 1985) indicated that CDM with this normalization 
was not a good match to the galaxy redshift data on scales of a few Mpc, in particular 
the galaxy pairwise velocities <7i 2 (l/j _1 Mpc) 340 km s (D avis & Peebles 1983, hereafter 
DP83). CDM is a better fit to this data for a 8 « 2.5. 

As more data became available on larger scales, however, it became increasingly clear 
that this data favored a higher normalization (see e.g. the reviews by Davis et al. 1992, 
Ostriker 1993). Thus, the problem with CDM is — at least — that the spectrum shape is 
not quite right. It is of course also possible that the problem is much more serious: that the 
dark matter in the universe is not mostly cold plus a little baryonic dark matter (presumably 
relatively cool ionized gas and MACHOs). But since there are at least two well-motivated 
candidate cold dark matter particles (axions and lightest superpartner particles; see e.g. 
Primack, Seckel, & Sadoulet 1988, Griest 1995), and since it is remarkable that the simple a 
priori CDM model came so close to fitting all the data, we will assume that the dark matter 
is indeed mostly cold, but that there is some aspect of the theory that needs to be modified 
slightly. Many variants of CDM were studied by Holtzman (1989) and Holtzman & Primack 
(1993) and others (for a review see Liddle & Lyth 1993). Of all the variants of CDM with 
= 1, the best fit to the data appears to be a model where, along with mostly cold dark 
matter, there is also some hot dark matter (light neutrinos). 

We call this theory Cold + Hot Dark Matter (CHDM). The first detailed simulations of 
CHDM showed that it indeed appears to be a good fit to all available present-epoch data, 
including both large and small-scale galaxy velocities (Klypin et al. 1993, hereafter KHPR; 
cf. Davis, Summers, & Schlegel 1992; Jing et al. 1994; these papers give references to earlier 
papers on CHDM). 

In the present paper we describe a series of higher-resolution simulations that we have 
performed in order to understand better the differences between CDM and CHDM, in par- 
ticular, the details of the distribution and velocities of dark matter halos. In this paper, we 
describe the simulations including the evolution of the power spectra (§3), and our simplest 
schemes for identifying "galaxies" in these simulations (§4). We then show that the corre- 
sponding CHDM galaxy power spectra are in good agreement with recent observations of 
both CfA2 optical and IRAS galaxy redshifts on both large and small scales (§5); and we 
compare correlation functions of galaxies and cold dark matter particles in real space, and of 
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galaxies compared to observational data in redshift space (§6). Finally, we discuss velocity 
dispersion and velocity bias in our simulations (§7), and summarize our conclusions (§8). 

In a series of related papers, we have further analyzed these same simulations in a variety 
of ways. In Yepes et al. (1994), the results were compared to observations of the galaxy angular 
correlation function; this paper also includes a discussion of the time evolution of correlations 
in the simulations. In Nolthenius, Klypin, & Primack (1994, NKP94), we presented the 
highlights of our analysis of the velocity structure of the galaxy groups in these simulations, 
compared to the CfAl data. We merged the CfAl galaxies that we could not have resolved, 
determined the corresponding Schechter luminosity function, assigned Schechter-distributed 
luminosities in rank order to the halos in our simulation volume, and constructed magnitude- 
limited redshift-space catalogs that we analyzed in parallel with the merged CfAl catalog. 
Our conclusion was that CDM groups have systematically higher median velocity dispersion 
than the CfAl data, as a function of either the redshift- link parameter or the fraction of 
galaxies in groups. Visualizations of these simulations — including a videotape — show that 
the CHDM simulations are much more filamentary than the CDM simulations (Brodbeck et 
al. 1995, hereafter BHNPK). Nolthenius, Klypin, & Primack (1995, NKP95) presents a more 
detailed group analysis, including a number of alternative methods of assigning luminosities 
and identifying galaxies, and several schemes for breaking up overmerged halos. These papers 
show that CHDM with fi c /fij, /fi& = 0.6/0.3/0.1 is a much better match to this data than 
CDM, but the more detailed analysis of NKP95 suggests that the fit would be somewhat 
better with a little lower fij,. This also appears to be indicated by the analysis of the Void 
Probability Function (VPF) for these same simulations by Ghigna et al. (1995), which finds 
that our CDM simulations are a better match to the Perseus-Pisces Survey data than our 
CHDM simulations. 

Several papers (Kauffmann & Chariot 1994, Mo & Miralda-Escude 1994, and Ma & 
Bertschinger 1994) pointed out that CHDM with fij, = 0.3 does not produce enough dark 
matter halos to account for the data on damped Lyman a systems at redshifts z > 3. We 
(Klypin et al. 1995, KBHP) agree with this (although our calculations disagree in some 
details with those of the above papers), but we show that if fij, is lowered a little — to 
fijz ~ 0.2, corresponding to ~ 5 eV of neutrino mass for h = 0.5 — there are enough halos 
with characteristic velocities < 100 km s _1 to correspond to the observations, as long as 
the number of these damped Lya systems decreases above about z > 3.5. Recent neutrino 
experimental data suggests that ~ 5 eV of neutrino mass is shared about equally among 
two species of neutrinos, v m u and v Tl and linear calculations suggest that this will be an 
even better fit to the available data than CHDM with one 5 eV neutrino (Primack et al. 
1995). We are presently analyzing both dissipationless and hydrodynamical simulations with 
two neutrinos, using both the NKP analysis and several new statistics — including several 
varieties of shape statistics (moment statistics: Hellinger et al. 1995; filament statistics: Dave 
et al. 1994, 1995; etc.) and also an improved version of the Davis-Peebles (1983) pairwise 
velocity statistic (Somerville, Davis, Nolthenius, & Primack, 1995) — which we have been 
developing using the suite of simulations described in the present paper as a test-bed. 

2. COSMOLOGICAL MODELS 

We describe simulations of two cosmological models in this paper: Cold + Hot Dark 
Matter (CHDM) and Cold Dark Matter (CDM). Both models assume that the Universe has 

critical density (O = 1), that the Hubble parameter is h = 0.5 (H = 100/i km s" 1 Mpc -1 ), 
and that the cosmological constant is zero. The CHDM model has the following parameters. 
The mass density of the Universe in the form of light neutrinos is Q v = 0.30, and the density 
of baryons is = 0.10. We use the analytical approximations eq. (1) of KHPR for the "cold" 
and "hot" spectra for the CHDMi simulation, and the spectra given in KHPR Notes added 
in Proof (eq.Al-2) for the CHDM2 simulation. The amplitude of fluctuations is normalized 
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so that our realizations are drawn from an ensemble producing the quadrupole in the angular 
fluctuations of the cosmic microwave background at the Q ps -norm = 17 /iK level measured 
by COBE (Smoot et al. 1992). In linear theory, this corresponds to rms fluctuations of mass 
in a sphere of 8/i -1 Mpc radius <7g = 0.667 and the rms velocity of the "cold" dark matter 
relative to the rest frame a v = 750 km 



s 



Our CHDM model perhaps assumes an unrealistically large fi&. A lower value = 0.05 
is better consistent with standard Big Bang Nucleosynthesis (Walker et al. 1991), although 
the large X-ray luminosity of rich clusters (oc fijj) suggests higher (e.g. White et al. 1993) 
if O = 1. Recent BBN studies also suggest a higher upper limit < 0.1 for h = 0.5 (Copi, 
Schramm, & Turner 1995). The power spectrum of fluctuations in the CHDM model increases 
by approximately on 20% scales of < 10/i _1 Mpc if we change from 0.1 to 0.05 and keep 
the same normalization on very large scales. Thus the results presented in this paper are also 
approximately valid for a CHDM model with = 0.05 and the amplitude of the quadrupole 

Qps— norm — 14 fili. . 

Another possible interpretation is a model with = 0.05, Q ps -norm = 19 /J.K as indi- 
cated by the latest COBE analyses (Gorski et al. 1994, Bunn, Scott, & White 1995), with 
a 10% contribution of gravity waves to the AT/T fluctuations on COBE scales and a 20% 
further reduction of power on scales of ~ 20/i _1 Mpc due to a small amount of tilt — i.e., a 
primordial power spectrum P{k) oc k n with n 0.95 instead of the Zel'dovich value n = 1. 
This amount of gravity waves and tilt is consistent with expectations for simple inflationary 
models with inflaton potential V{cf>) = \cf> a with a = 4 (see e.g. Turner 1993). Our CHDM 
normalization is in any case within the COBE errors. 

We use the spectrum given by Bardeen et al. (1986, BBKS) for the CDM model. Two 
normalizations are used for CDM: a "biased" model — which we call CDM1.5 — with linear 
bias factor b = a^ 1 = 1.5 or <7g = 0.667 (which leads in linear theory to a v = 660 km s _1 
and Qps-norm = 8.5 /iK), and an "unbiased" CDM model — called CDM1 — with <7g = 1, 
o v = 990 km s _1 , Q ps -norm = 12.8 /iK. This CDM1 model normalization was compatible 
with the lowest normalization quoted from the first-year COBE measurements (Smoot et al. 

1992) , but in view of the more recent COBE results it would perhaps be more reasonable 
to interpret this model as corresponding to about 35% contribution of gravity waves to the 
COBE amplitude but with no tilt in the spectrum. This is not what is expected in most 
cosmic inflation models, but "natural inflation" models (e.g. Adams et al. 1993, Turner 

1993) do produce gravity waves without tilt. 

It is perhaps worth noting that various CDM spectra are not entirely equivalent (cf. also 
Liddle & Lyth 1993). For example, Efstathiou, Bond, & White (1992), using the spectrum of 
Bond & Efstathiou (1984), quote Q ps -norm = 15.3 /iK as implying erg = 1.10 for CDM; thus, 
Qps-norm = 17 /iK would correspond to erg = 1.22. With the BBKS spectrum, Q ps -norm = 
17 /iK corresponds to erg = 1.33, or b = 0.75, for CDM. 



3. N-BODY SIMULATIONS 

Numerical simulations were done using a standard Particle-Mesh (PM) code (Hockney 
& Eastwood 1981). The equations we actually solve are given by Kates, Kotok, & Klypin 
(1991). The code is very fast and it is especially suitable for analysis of statistical properties 
of cosmological objects (galaxies, groups, and clusters). Two conditions should be satisfied 
at the same time for this kind of analysis: the volume of the computational box should 
be representative for the objects under investigation, and the resolution should be good 
enough to allow the identification of the objects and to allow us to estimate a few basic 
parameters of the objects (e.g., the mass distribution and the velocity relative to the rest 
frame) in the simulation. There are numerical methods which can provide better resolution 
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of the gravitational forces (e.g., P 3 M or tree codes; for a review, see Sellwood 1987), but the 
number of objects and the mass resolution — and thus the resolution of the initial spectrum 
of fluctuations — is worse compared to what the PM code provides. Those methods are more 
suitable, and at present are more often used, to study separate objects or a few objects (e.g., 
the collision of galaxies, the formation of a galaxy cluster). We have had experience with such 
codes, and also with hydrodynamics codes, but for our present purpose PM code is superior. 
This is especially true for CHDM, for which the high velocities of the hot dark matter require 
that we have many hot particles to adequately sample the six-dimensional phase space. 

This paper discusses three runs: two for the CHDM model (CHDMi and CHDM2) and 
one for the CDM model, but we analyze results from the CDM run at two different expansion 
factors, corresponding to linear bias parameters b = 1.5 (CDM1.5) and 6=1 (CDM1). All 
simulations were done using a 512 3 force mesh. The CHDMi and CDM runs were started 
with fluctuations generated from the same set of random numbers. 

Each simulation has 256 3 "cold" particles. The CHDM simulation has two additional 
sets of 256 3 particles to represent the "hot" neutrinos. We use the same prescription to 
simulate random thermal velocities of the "hot" particles as KHPR. "Hot" particles are 
generated in pairs, particles of each pair having random "thermal" velocities of exactly equal 
magnitude but pointing in opposite directions. The directions of these "thermal" velocities are 
random. The magnitudes of the velocities are drawn from relativistic Fermi-Dirac statistics 
(see KHPR). The CHDM spectra we use were calculated treating cold dark matter, hot dark 
matter (neutrinos of 7 eV mass), and baryonic matter separately. However, as usual in CDM 
simulations, when we switch from linear to nonlinear evolution, the cold particles represent 
both cold and baryonic dark matter. The particles have different relative masses: each "cold" 
particle has a relative mass 0.7 and each "hot" particle has relative mass 0.3/2 = 0.15. 

The size of the computational box for both the CHDM and the CDM simulations is 
100 Mpc (i.e., 50/i _1 Mpc for h = 0.5). The smallest resolved comoving scale in these simula- 
tions is 97.6/i _1 kpc and the mass of a "cold" particle was 1.45 X 10 9 h~ x Mq for the CHDM 
simulation and 2.07 X 10 9 h -1 Mq for the CDM simulation. Both CHDM simulations were 
started at redshift 15 and were run to redshift zero with a constant step Aa in the expansion 
parameter a. The CDM simulation was started at z = 18 for the biased model and z = 27.5 
for the unbiased model. The step for CHDMi and CDM simulations was Aa = 0.01; for the 
CHDM2 simulation it was a factor of two smaller, Aa = 0.005. 

After running the CHDMi simulation, we found that there were two mistakes in our 
initial conditions (see KHPR, Note Added in Proof): the fitting formula for the cold spectrum 
was too small and the velocities were too large, both by about 20% on small scales. However, 
these effects are in phase and largely cancel. We have run another CHDM2 (rev) simulation 
with both errors corrected, and found that the power and velocity differences between this and 
CHDM2 (old) on small scales declined to 5% by z = 7 and remained at this level. All results 
discussed here are from CHDM2 (rev), but the agreement with CHDM2 (old) is well within 
the la error bars; thus the CHDMi simulation should also be reliable. (This is discussed 
further below in §5.) 

Initial positions and velocities of particles were set using the Zel'dovich (1970) approx- 
imation. The displacement vector was simulated directly. Phases of fluctuations were ex- 
actly the same for "hot" and "cold" particles. When generating velocities of "hot" particles, 
the "thermal" component, as described above, was added to the velocity produced by the 
Zel'dovich approximation. In the case of the CHDM2 (rev) simulation, (referred to hereafter 
as CHDM2) the Zel'dovich approximation was slightly modified to take into account the fact 
that the growth rate of fluctuations is different for different wavelengths (Ma 1993). The 
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correction was made as follows. If is the spectrum of the velocity potential at the ini- 
tial moment and a is a normalization constant, then the relations between the Lagrangian 
(unperturbed) position q, Eulerian position x, and peculiar velocity v are given by 

x = q — aa(t) k$yt exp(— zkq), v = adx/dt = — aaa fk^^k exp(— zkq), (2) 

where fk = din S co \d(k) / din a is given by eq.(A2) of KHPR (for z = 15 fk changes from 1. 
at small k to 0.805 around k = 0.3Mpc _1 ). Wh en setting up initial conditions for CDM and 
CHDMi we used f k = 1. 

When we completed the CHDMi and the CDM simulations we realized that there was 
a statistical fluke of probability ~ 10% affecting the longest waves in the simulations: the 
amplitude of the waves was a factor of 1.3-1.4 larger (so the power was about a factor of 2 
larger) than that expected for the ensemble average. Nevertheless, the waves (A = 100 Mpc) 
are still in the linear regime even at the end of the simulations and it is relatively easy to 
make corrections to the correlation functions and the power spectrum. Moreover, as we argue 
below, this small additional power could be considered as a compensation for the finite size 
of the simulations. It brings the level of fluctuations inside our 100 Mpc box close to that 
expected for a much larger box. The CHDM2 simulation has power typical for a box size of 
100 Mpc. Differences between the CHDM simulations — which are readily apparent in our 
BHNPK visualizations — thus indicate to what extent various statistics that we use in this 
and other papers are affected by cosmic variance. 

The CHDM2 simulation does not have that excess power on long waves and, as a result, it 
does not have a structure comparable with the Great Wall, which to some extent dominates 
the CfA2 (m < 15.5) catalog. Nevertheless, not all statistics are strongly affected by the 
presence of the large structure. We found that the pair-wise velocities and fraction grouped 
were affected the most (because of a larger fraction of galaxies being in galaxy groups with 
higher velocity dispersion). The power spectra and group velocities were less sensitive to the 
effect (NKP95). 

4. BIASING SCHEMES AND GALAXY FINDING ALGORITHMS 

We use the following procedures to identify "galaxies" in our simulations in this paper. 
The density field is produced on the 512 3 simulation mesh. We find maxima of the total 
("cold" plus "hot") density. Then all local maxima above the density threshold /? t hr > 30/9 cr it 
are found. This list played a role of a pool from which more dense objects are chosen for data 
analysis. The overdensity threshold corresponds to a mass of 7.76 X 10 9 h~ x Mq in a cell. The 
number of maxima depends on the simulation considered and the moment of time at which we 
consider it. For CHDMi and for the "unbiased" CDM simulations at the final moment there 
are 29151 and 37164 dark halos respectively. Usually we use more massive objects selected 
from the lists of dark halos, with the mass limits 1.27 X 10 10 h~ x M Q and 2.5 X 10 10 h~ x M Q 
being typical choices. Because the mass function of dark halos n(M) decreases with mass M 
approximately as n(M) oc M~ 2 (implying that the cumulative number N(> M) oc M _1 ), 
doubling the mass limit reduces the number of "galaxies" by about a factor of two. (For more 
details, and results using several alternative galaxy identification and luminosity assignment 
schemes, see NKP95.) 

We use two estimates of the dark halo mass: the mass m\ in the cell of the density 
maximum, or the mass in the cube of 3 X 3 X 3 cells centered on the maximum. Effectively 
these definitions correspond to (i.e. have the same volume as) spherical "galaxies" of radii 
61/i _1 kpc and 182/i _1 kpc, respectively . In order to find the velocities of these dark halos, 
the velocity field was constructed on three 512 3 grids (one for each of the vector components) 
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using the Cloud- In-Cell technique. The velocity at the position of the density maximum was 
assigned to the dark halo. 

Figure 1 and Figure 2 present examples of two groups of galaxies identified in the 
CHDM2 simulation. Groups were found as local maxima of "cold" mass within sphere of 
0.75/i -1 Mpc. Positions of all cold particles inside a 3/i -1 Mpc cube are shown together with 
positions of dark halos with central overdensity larger than Sp/p > 150. Each circle on the 
plots corresponds to a dark halo, with the area of the circle being proportional to the mass 
of the dark halo. The group in Figure 1 is quite massive: M = 7.5 X 10 13 h~ x Mq . There 
are about 66000 cold particles in the frame. The group has nine "galaxies" with Sp/p > 150. 
("Galaxies" with X > 89.5 Mpc do not belong to the group). The cold particle rms 3D 

velocity for that group is 644 km s -1 . Figure 2 shows a smaller group: the cold particle rms 
3D velocity is 310 km s" 1 , and M = 1.85 x 10 13 // '.\/. . The group has four members of 
about equal mass. Five "galaxies" on the upper part of the plot and one small dark halo at 
the bottom are not the members of the group. 

5. POWER SPECTRA AND CORRELATION FUNCTIONS 

Initial spectra for the simulations are shown in Figure 3 . The spectra are given at 
z = 15. The CDM simulation was started earlier (z = 18), but its initial spectrum was 
scaled up by linear theory to the moment z = 15. We present spectra up to the wavenumber 
corresponding to the Nyquist frequency of "particles: k = (2tv / L) X (iV part /2) = 16h Mpc -1 , 
where L = 50/i -1 Mpc and iV part = 256. (Recall that we have 256 3 cold particles in each 
simulation.) Because of the small number of harmonics with long waves, statistical fluctua- 
tions are apparent relative to theoretical spectra at k < lh Mpc -1 . Nevertheless, the overall 
accuracy of normalization is better than 5%. 

At k > lh Mpc -1 there was 10-20 percent difference in the power spectra between the 
CHDMi and CHDM2 models. The difference originates in different approximations used for 
the simulation. The CHDMi simulation was done for the spectrum given by eq.(l) of KHPR. 
The spectrum actually was the power spectrum of baryons, not cold particles. As one might 
expect, the difference between the spectra is small because baryons had enough time from the 
recombination to z = 15 to fall in the potential well defined by the dark matter, but 10-20 
percent on the smallest scales is not negligible nowadays. Another effect on the same level 
basically compensates the difference in the spectra. When setting initial conditions for the 
CHDMi simulation, we forgot to take into account that the growth rates of fluctuations are 
slightly smaller as compared with pure CDM case on scales smaller than the effective Jeans 
length in neutrinos (see eq.(2)). Thus velocities were set 10-20 higher than they should be. 
Both effects - smaller densities and larger velocities - have the same phases and amplitudes, 
but opposite signs. So after few expansions the system (still in linear regime) arrives at almost 
the same configuration as if it started from correct initial conditions. 

The CHDM2 simulation was done for the spectrum (Al-2) of KHPR, which does not 
have the problems just described for CHDMi . In order to check the effect of the difference 
in spectra and velocities, two tests were done. One was discussed in detail by KHPR in 
Note added at proof. Two models with the same realization of random numbers were run up 
to z = 1.48, one model having 10% higher initial density amplitude and 20% lower initial 
velocities than the other. At the final moment positions and velocities of particles in the two 
simulations were compared one-by-one. The rms difference in particle positions was 0.11 cell 
units and the rms velocity difference was 1.6%, i.e. truly negligible. Another test involved 
comparing the CHDM2 simulation with one run with the same initial random numbers, but 
with the initial conditions set exactly as for the CHDMi simulation (that is, eq. (1) of KHPR 
for the power spectrum and no corrections for the growth rates). By redshift z = 7 the 
difference in power spectra of the two simulations became smaller than 5%. The difference 
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between the runs at the end of the simulations (z = 0) was also small. For example, the 
number of dark halos with mass larger than M > 1.8 X 10 10 // '.\/. was 29662 in the 
"uncorrected" simulation and 29795 for the CHDM2 simulation. Even the numbers of the 
most massive, rare dark halos were quite close: 88 and 81 for M > 5 X 10 12 h~ x Mq . (All 
masses are estimated for 3 3 cells centered on density maximum.) 

In Figure 4 , we present the evolution of the power spectrum in the CHDMi model. The 
three solid curves show the power spectra of the total density ("cold" plus "hot") in the 
simulation at z = 15 (the lowest curve: initial moment), at z = 3.70 (the middle curve), 
and at the present moment (z = 0, the top curve). The dot-dashed curves present the linear 
theory predictions. Effects of nonlinearity are obviously seen in the spectrum at z = for 
k > 0.2 h Mpc - . First, there is a long tail with slope n ~ —1.2 from k = 0.2 h Mpc - to 
k = 2 h Mpc -1 , which is responsible for the power-law behavior of the correlation function. 
Second, at higher wavenumbers (k > 2 h Mpc -1 ) the power spectrum has a tail P(k) oc A; -2 ' 5 , 
which agrees with the BBGKY similarity solution for clustering in statistical equilibrium 
(Peebles 1980, §73). (The BBGKY solution predicts the slope —6/(5 + n), where n is the 
initial slope of the power spectrum, which for the CHDM model in this range of wavenumbers 
is n = —2.6.) As we mentioned above, due to a statistical fluke the first harmonic in the 

CHDMi and CDM simulations (k = 0.125 h Mpc -1 ) is above the linear prediction just at 
the beginning of simulations. But it stays in the linear regime even until z = 0. In Figure 5 
we compare the power spectra in the CHDMi and the CDM simulations. The spectrum of 
the CHDMi and that of the biased CDM (erg = 0.66) simulations are quite close in the range 

k = (0.1 - 1) h Mpc -1 , but their long waves are very different: the biased CDM spectrum is 
higher by more than a factor of three. Also the slope of the power spectrum in the nonlinear 
regime k = (0.2 — 2.0) h Mpc - is shallower for the CDM models (n ~ — 1.) as compared 
with the CHDMi spectrum, which results in a steeper correlation function at small radii for 
CHDMi. 

The power spectra obtained in the simulations are compared with observational data 
in Figure 6 and Figure 7 . Because observational spectra were measured in redshift space, 
we simulated the distribution of "galaxies" in redshift space s by placing an "observer" in a 
corner of the simulation box and by displacing dark halos along the line of sight in accordance 
with their velocity v and distance r from the observer: s = r(l + (vr) / ' Hr 2 ). Dark halos with 
central overdensity more than 100 were selected as "galaxies". Redshift distortions produce 
an effect that is analogous to Malmquist bias: some of the halos close to boundaries of the 
simulation box and moving away from the observer can disappear from the box, reducing the 
mean density at the boundaries. (Note that we cannot apply periodic boundary conditions for 
those receding halos). In order to avoid this Malmquist-like bias, we periodically replicated 
all halos before going to redshift space, but then only retained the halos inside the box. 
The mean number density of "galaxies" in real and redshift space was almost the same: 
0.0525/i 3 Mpc -3 in redshift space and 0.0508/i 3 Mpc -3 in real space. The density field of the 
halos was constructed using the Cloud- In-Cell technique on our 512 3 mesh and was smoothed 
using a Gaussian filter with 1 cell width to reduce the shot noise. This did not remove all the 
noise: there was apparent flattening of the spectrum at k > 2k Mpc -1 . A constant P = 12 
was subtracted from the power spectra to approximately remove the shot noise. In addition, 
the power in the first k bin (longest waves in the simulations) in the CHDMi simulation 
was reduced by two to take into account the excess of power in the initial spectrum for this 
simulation. Besides that difference in the first bin, both CHDMi and CHDM2 simulations 
produce remarkably similar power spectra. 

Optically selected CfA and IR-selected IRAS galaxies have different bias relative to each 
other, and thus sample the underlying power spectrum of the dark matter differently. So we 
clearly cannot compare both these data sets the same way to dark matter particles or halos. 
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Instead, we adopt what we think are simple and resonable prescriptions. When estimating 
the CHDM power spectrum (solid curve: CHDMi, dotted curve: CHDM2) to compare to 
that of the CfA2 galaxies (Vogeley et al. 1992) presented in Figure 6 , each dark halo in the 
numerical simulations was weighted by its mass (defined as mass within 3 3 cells centered on 
the peak of density). This weighting procedure was designed to compensate for the lack of 
resolution in the central parts of clusters, which leads "overmerging" at the centers of clusters 
— a very massive dark halo representing in reality several "galaxies," not one. To the extent 
that CfA-selected galaxies have the same mass-to-light ratio, by mass-weighting the halos 
we are mimicking this in the solid curve. It is less clear how to mimic the IRAS galaxies, 
presented in Figure 7 (dots — results of Fisher et al. 1992, dashed broken curve — results of 
Feldman et al. 1993). The solid and dotted curves attempt to mimic IRAS galaxies by equally 
weighting all dark halos with density greater than 100/? c . This prescription underestimates 
the central regions of clusters and large groups, just as IRAS galaxies do. 

Because of the quite limited size of our computational box, we cannot confront the results 
of our simulations with observational data on the power spectrum at small k. Nevertheless, it 
seems quite reasonable to suggest that on long waves the redshift distortions do not depend on 
the wavelength (Kaiser 1987). As we discussed earlier, the longest waves in our simulations 
still grow in the linear regime. Thus we can use the power spectrum given by the linear 
theory to make a prediction for the spectrum of galaxies in redshift space. The normalization 
constant for the power spectrum in redshift space was found by matching powers at the 
longest waves in the simulations. The power spectra estimated in this way are shown as 
dot-dashed curves in Figure 6 and Figure 7 . The linear spectrum was multiplied by the 
factor 2.75 in the case of CfA2 galaxies and by the factor 1.47 for the case of IRAS galaxies. 
If we assume that the relation between the redshift P s and real-space P r spectra is given by 
(Kaiser 1987) 

P s (k) = P r (k) ■ b 2 (l + 2/36 + 1/56 2 ), (3) 

where b is the biasing parameter, then biasing parameters for CfA2 and IRAS galaxies are 
6 c fa = 1-3 and &iras = 0.85, which are close to other estimates of the biasing parameters 
for those types of galaxies. Unfortunately, these estimates are of uncertain accuracy both 
because of uncertainties in the matching procedure and because it is not clear whether the 
redshift correction is truly scale independent at k around 0.1-0.2 where we adjust the spectra. 
Nevertheless, we think it is quite plausible that these CHDM spectra agree reasonably well 
with power spectra determined from both CfA2 and IRAS galaxy catalogs. 

6. CORRELATIONS 

In Figure 8 , the correlation functions £ of different populations in the CHDMi model are 
shown at the present moment z = 0. The lower solid curve presents the correlation function 
of the "cold" dark matter in the simulation. In the range R = (0.5 — 5)/i _1 Mpc the correlation 
function follows the power law (,(R) = (i?/4.0/i -1 Mpc) -1 ' 8 indicated by the lower dot-dashed 
line. At larger radii the correlation function goes above the power law, but reasonably close to 
the correlation function predicted by the linear theory for an infinite box (the dashed curve). 
Actually on scales around 10/i _1 Mpc, the level of fluctuations in the simulation is closer to 
that expected for much larger box-size. The correlation function crosses zero at 25/i _1 Mpc 
due to the finite size of the computational box. 

The correlation function of dark halos in the CHDM model is larger than the correlation 
function of the dark matter. The three upper solid curves on the plot show the correlation 
functions of dark halos satisfying three different mass limits. From the bottom to the top, the 
curves are for the following mass limits and numbers of "galaxies": m\ > 1.28 X 10 10 h~ x Mq , 
N = 14606; m x > 2.56 x 10 10 Ii^Mq , N = 6172; and m } > 7.67 x 10 10 Ii^Mq , N = 1410. 
Here m\ indicates that the mass corresponds to the mass in one cell at the density maximum. 
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The mass of a halo typically increases by 7-10 times when one goes from m\ to the mass 
defined in 3 3 cells. When we estimated the correlation function, a dark halo was weighted 
proportionally to its mass found in the 3 3 cells volume centered on the density maximum. The 
reason we use this weighting of dark halos is to compensate for the effects of "overmerging" 
of dark halos in central regions of groups and clusters. 

The growth of correlations when the mass limit increases is apparent, but it is not large. 
At radii smaller than 5/i _1 Mpc the correlation function is approximated by the power law 
(,(R) = (-R/6.7/i -1 Mpc) -1 ' 8 . This corresponds to the bias parameter b = \J £ ga i/ Cdm ~ 1.6 
in a reasonable agreement with the linear bias &ii nea r = l/cs = 1-5. The correlation length 
increases by about 1 Mpc when the number of "galaxies" drops by a factor of ten because we 
have selected a higher mass limit. 

Any type of N-body simulations do not include waves longer than the size of the initial 
cloud of particles. We hope that those long waves do not change our results too much. 
Different statistics are sensitive in different degree to those waves. By comparing results 
from simulations with different sizes or variations from one realization to another, one can 
get an estimate of the effect. Sometimes it is possible to make corrections, which take into 
account effects of the waves. Because the waves of the size of the computational box and 
larger are still in the linear regime, it is possible to make corrections for those waves and for 
the statistical fluctuations of the largest waves in the box. The relation between the power 
spectrum and the correlation function could be split into two parts: the contribution of waves 
in the simulation box and the contribution of waves longer then the length of the box. We 
have 

, 1 [ k ™* ,s\n(kr) „ 1 f k bo* sin(&r) „ , , 

£(r) = — / k 2 P(k) ) J dk + — k 2 P(k) ; J dk, (4) 

^ } 2vr 2 A box [ ) kr ^ 2vr 2 J kr ' [ ) 

where &box = 2tt/L } and & max = 256&box are the smallest and the largest wave numbers in the 
simulations (L = 50/i _1 Mpc is the box size). The first integral in this expression is estimated 
by the usual procedure of counting pairs of objects at given separation. The second part is 
estimated analytically and added to the first one. In order to make the correction due to 
the fluctuation in the initial spectrum, we estimate the additional power due to the fluke 
in the initial spectrum, extrapolate it by linear theory to the present moment of time, and 
subtract the contribution due to the power from the correlation function. This procedure is 
applied to radii smaller than 20/i _1 Mpc. At larger radii we use the linear theory. A test for 
the procedure is provided by comparing predictions of the linear theory with the corrected 
correlation function at radii (10 — 20)/i _1 Mpc, where the correlation function is small and 
one can hope that the linear theory predictions are still reasonably accurate. 

In Figure 9 , the correlation function of the "cold" dark matter in the CHDMi simulation 
is shown as the lower solid curve. The dashed curve presents the linear theory prediction. The 
agreement with the corrected correlation function in the simulation looks quite reasonable. 
The same procedure was applied to the correlation function of dark halos with m\ > 1.28 X 
10 10 h~ x Mq . In this case the power spectrum predicted by the linear theory was multiplied 
by a factor 2.3 to account for the biasing of dark halos observed in Figure 8 . The dot-dashed 
lines on the plot are the same power-laws as in Figure 8 . The CHDM model predicts that in 
real space the correlation function of both the dark matter and the galaxies is a power law 
on scales smaller than 20/i _1 Mpc with the slope 7 = —1.8. The correlation length of the 
dark halos is r ga i = 6.7/i _1 Mpc. The correlation function crosses zero at ro = 50/i _1 Mpc and 
becomes negative at larger radii. 

For comparison, Figure 10 presents the correlation function in real space for the unbiased 
CDM simulation. The correlation function of the dark matter particles is not a power law, as 
was well known (e.g., Suto & Suginohara 1991, Bahcall, Cen, & Gramann 1993). Nevertheless 
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the correlation function of dark halos does show power-law behavior (the dot-dashed line) 
with —1.8 slope and the correlation length r ga i = 5.7/i _1 Mpc. The same mass threshold 
(mi > 1.28 X 10 10 /i _1 M ; N = 22275) and weighting scheme was used as for the CHDM 
halos. No corrections were applied in this case. 

In order to estimate the correlation function in redshift space two "observers" were 
placed in opposite corners of the computational cube. Then objects of interest (dark halos or 
a small random fraction of the dark matter particles) were 27 times periodically replicated in 
space in such a way that the original box was surrounded by 26 replicated boxes. For each 
observer the objects were displaced along the line of sight in accordance with their peculiar 
velocities. Then neighbors at a given distance were counted for all objects inside the original 
box. Results of the two observers were averaged. 

Figure 11 shows the correlation function for the dark matter and dark halos in the 
CHDMi simulation. Note that at small radii the difference between halos and dark matter 
becomes larger, which indicates that the peculiar velocities of the halos are smaller than 
those of the dark matter, i.e., they exhibit a global velocity bias. (We saw this directly in 
KHPR; see also NKP95.) These peculiar velocities have two opposite effects on the correlation 
function: at small radii it became smaller, and at large radii it was increased. This agrees 
with theoretical expectations. Peculiar velocities on small scales produce a "finger of god" 
effect, which results in less compact clusters and groups. On large scales, where fluctuations 
are still in the linear regime, peculiar velocities displace particles torward collapsing regions, 
thus increasing apparent clustering. 

In redshift space corrections to the correlation function due to waves longer than our 
simulation box or due to statistical fluctuations of the longest waves in the simulations can 
be done in the same way as for £(r) in real space. In order to take into account the effect 
of peculiar velocities at large radii on we apply the Kaiser (1987) correction to the 

linear power spectrum given by (3) (note that the biasing parameter is 6 = 1 for the dark 
matter — it is unity even for "biased" models). For dark halos in the CHDM model the 
parameter is 6 = 1.5 in accordance with the ratio of correlation functions of "galaxies" and 
the dark matter in real space. For the unbiased CDM model we take 6=1. The corrected 
correlation functions for the CHDMi model are shown in Figure 12 as the solid curves. The 
dashed curve presents the linear theory prediction amplified by the velocity effects. Triangles 
show the estimates of the correlation function for the CfA Redshift Survey (Vogeley et al. 
1992). The agreement between the observational data and the model is reasonably good on 
all scales. There is a slight indication that the model predicts larger but taking into 

account all uncertainties involved on both the theoretical and the observational sides, the 
difference hardly is significant. 

The correlation function in redshift space for the unbiased CDM model is shown in 
Figure 13 (with no corrections) and in Figure 14 (corrected for the box size and Kaiser's 
correction for 6=1, with no corrections for the additional fluctuations in the box). The 
dashed curve in the last plot is for linear theory in redshift space. It really does not match 
well with from the simulation. Many different parameters were tested and this is the best 
match we can get. The triangles again show results of Vogeley et al. (1992). The correlation 
functions of the dark matter and the dark halos ("galaxies") match very well at scales smaller 
than s < 3/i _1 Mpc, but £ ga i is about 1.5 times higher than ^ m for s < 3/i _1 Mpc. This 
is probably the reason why the linear theory prediction does not match the results of the 
simulations: larger 6 would increase £iinear(s), but unfortunately that would contradict the 
very tight match of ^ m and £ ga i in real space at this radius. In any case, the unbiased CDM 
model predicts too small £ ga i(s) at s < 5/i _1 Mpc and it fails to get the observed level of 
the correlation function at scales around 30/i _1 Mpc. As a matter of fact, no matter what 
amplitude one assumes or what correction one applies, the CDM model cannot produce the 
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observed correlation function at ~ 30h 1 Mpc because at that radius the CDM correlation 
function crosses zero. 

7. VELOCITIES 

Velocities provide very important but somewhat controversial information. One of the 
most important issues is the pair-wise velocities on small scales. For a long time the results 
of DP83, who found that the rms galaxy pairwise velocity (which they designated as a(r p )) 
is <7i2(l/i -1 Mpc) = 340 km s— 1, have been a stumbling block for cosmological models. One 
should mimic all the steps of DP83 in numerical simulations in order to obtain directly 
comparable results. In KHPR we only partly mimicked the DP83 analysis; in particular, 
their infall model was not simulated. We showed that the CHDM model predicts smaller 
pair-wise velocities than the CDM model, and that although the dark matter particles have 
higher velocities, the pairwise velocities of dark halos are consistent with DP83. In a similar 
calculation, Cen and Ostriker (1994) found somewhat higher pairwise velocities in simulations 
with the same CHDM parameters, and also found that the velocities depend on the local 
density. In fact, because the pairwise velocity statistic (T\2 is pair weighted and clusters not 
only have many close pairs but also high relative velocities, (T\2 is extremely sensitive to the 
presence of clusters in the volume covered by the redshift survey (da Costa, Vogeley, & Geller 
1994) or simulation; it is thus not a very robust statistic. Moreover, other authors who have 
recently tried to reproduce DP83 (T\2 results from the same CfAl data have been unable to 
do so (Mo, Jing, & Borner 1993; Zurek et al. 1994; Somerville et al. , 1995). It turns out that, 
due to a typographic mistake in the computer code used in the (T\2 calculation, the entire 
core of the Virgo cluster was inadvertently omitted; when it is included, <7i2(l/i -1 Mpc) 

600 km s _1 . The fact that IRAS galaxies give smaller u\2 (Fisher et al. 1994) doubtless 
reflects the fact that they are not found in clusters. 

We do not repeat the DP83 analysis in this paper, and we postpone detailed comparison 
of pair- wise velocities with observation for a later publication. (As mentioned in §1, we have 
compared another velocity-related statistic, the median group velocity from the CfAl data, 
with the present CDM and CHDM simulations in some detail in NKP94 and NKP95, and 
concluded that CHDM with fi c /fij, /fi& = 0.6/0.3/0.1 probably has velocities that are a little 
smaller than the data.) Instead we study here another important and controversial statistic: 
the velocity bias and real-space pair-wise velocities. 

Figure 15 presents one of the components of the real-space pair-wise velocity tensor - 
rms velocities along the line joining the pair of objects or cry . The top panel shows cry (R) for 
dark halos with central overdensity larger than 100. The mass inside the cell centered on the 
maximum is m\ > 2.5 X 10 10 h~ x Mq . The dot-dashed curves are for the CDM simulations 
(the higher one corresponds to b = 1.0, the lower one is for b = 1.5). The CDM b = 1.5 
results are close to CHDMi on scales around l/i _1 Mpc, but 6=1 CDM shows significantly 
larger pair-wise velocities. Note that the CHDM and CDM models have different trends. The 
CDM curves have maxima at 1.5 — 2/i _1 Mpc and at larger radii cry decreases. The CHDM 
curves do not show maxima in <7||(i?), but the slope of the curves is significantly smaller at 

R > 1.5/i _1 Mpc than at lower R. This difference is probably related to the difference in the 
large-scale power in the two models. 

The bottom panel shows results for the CHDM2 simulation. The three solid curves 
correspond to different mass thresholds: the larger the mass of both "galaxies," the larger 
their relative velocity. This trend can be explained if a large fraction of mass on these scales 
is associated with the dark halos. Then the more massive the halos that are selected, the 
deeper is the potential well and the larger are the relative velocities. Artificial "overmerging" 
due to the lack of either force or mass resolution would give another trend, with more massive 
halos moving slower. We do not see this trend in our simulations. A significant "overmerger" 
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problem appears in our case at overdensities larger than 1000. For comparison, we also show 
in Figure 15 the pair-wise velocity dispersion from KHPR simulations for the same CHDM 
model. The lack of long waves in the smaller box used by KHPR (L = 25/i _1 Mpc) resulted 
in smaller cry, but the agreement on scales less than l/i _1 Mpc is good. 

The tendency for more massive halos to move faster in pairs is not because all massive 
halos are moving faster; there is such tendency, but it is far too weak. The effect appears 
only in pairs of "galaxies." In Table 1 we present 3D rms velocities relative to the rest frame, 
a v , for different populations of objects. For the CDM simulation there is visible difference 
between a v of the dark matter and that of dark halos: "galaxies" are moving slower than 
dark matter particles in the CDM model. Nevertheless, the difference between galaxies with 
different masses is barely noticeable. The same is true for CHDM simulations, but in this 
case there is almost no difference between the rms velocities of cold dark matter particles and 
those of "galaxies." 

The ratio of rms velocities of dark halos to that of dark matter (both measured relative 
to the rest frame) can be called the global velocity bias: 



y\oh 



cr^ (galaxies) 

(5) 



cr^(dark matter) 



Table 1, fifth column, gives 6 U)g i b for our simulations. It is quite straightforward to add 
corrections to b v due to waves longer than the box size because those waves are still in the 
linear regime and they simply move our computational box as a whole. Column 7 and 8 give 
corrected a v and & U; gi b- The CDM model has global velocity bias & U; gi b ~ 0.8, while both 
CHDM simulations give numbers very close to unity: & U)g iob,CHDM ~ 1-0. 

In Table 1, columns 4 and 6, we present cry measured at R = l/i _1 Mpc and relate to it 
the pair-wise velocity bias: 

2 1 <7y (halo — halo) 

u ' pairs ^ ^ ^ erjj (dark matter — dark matter) ^ ^ 

The pair-wise velocity bias is significantly smaller than the global velocity bias and, for halos 
with the same mass, it does not depend on either the model or the large-scale power (com- 
pare CHDMi vs. CHDM2, and present paper vs. KHPR). Nevertheless, there is significant 
dependence of 6 U)Pa irs on the mass of dark halos. 

Thus, on small scales (r < 2/i _1 Mpc) the CDM model predicts pair- wise velocities larger 
than in CHDM. CDM normalized to COBE (b = 1.0) gives pair- wise velocities which may be 

too high (ss 600 km s 1 at l/i _1 Mpc) to be compatible with observations. On large scales 
r > 2/i _1 Mpc the CHDM model predicts cry which is between the CDM results for the model 
with the same biasing parameter b = 1.5 and the CDM with 6 = 1. It is interesting to note 
that for the CHDM cry is a rising function (at least up to 10/i _1 Mpc), which agrees at least 
qualitatively with observations (DP83, Mo et al. 1993). As we remarked above, in order to 
compare with observational results one needs to mimic the observations and the procedure 
for estimating the "observational" pair-wise velocities, which we are in the process of doing. 
Nonetheless, we can naively compare cry with observational results as other authors have 

done. According to Mo et al. (1993), the relative Id velocity in pairs with r p = l/i _1 Mpc 
is in the range 250 — 450 km s _1 . Here we take into account only optical samples (IRAS 
galaxies should be mimicked in a different way) and we do not consider the CfA2 sample with 
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Coma and other clusters for two reasons: i) The sample is heavily dominated by the Great 
Wall and is not a fair sample: the correlation function is far too big. ii) In our numerical 
simulations there are not such large clusters as Coma. Thus, it is not fair to compare the 
simulations, which do not have a cluster of that size, with a piece of the Universe, which 
does have a very large cluster. In the simulations cry is the range 350-450 km s _1 for the 
CHDM simulations, which probably indicates that CHDM is compatible with observations. 
The COBE normalized CDM with 6 = 1 600 km s 1 , which as we said is probably 

too large. 

Our results should be compared with results of other publications. Gelb, Gradwohl, & 
Frieman (1993) give cry = 1100 km s 1 at l/i _1 Mpc for dark matter particles in CDM with 

6=1, which is practically the same as we found in this paper (1160 km s _1 ). Because 
of their low resolution, Gelb, Gradwohl & Frieman do not give estimates for dark halos. 
Gelb & Bertschinger (1994) give a slightly larger number for the same CDM model: cry = 

1200 km s _1 . Their resolution was significantly better than that of Gelb, Gradwohl, & 
Frieman, and they also present estimates for pair-wise velocities of dark halos: 450-650 km 
s _1 at l/i _1 Mpc for "galaxies" with different circular velocities. Unfortunately, we cannot 
give reliable estimates for rotational velocities, but pair-wise velocities for our dark halos in 
CDM with 6=1 are in the same range as found by Gelb & Bertschinger. Nevertheless, 
Gelb & Bertschinger find the same trend as we do, for pair-wise velocities to be larger for 
more massive dark halos. There is also a reasonable agreement on the velocity bias in pairs: 
6 u , P airs(i? = l/i _1 Mpc) w 0.4 - 0.6. 

Comparison with KHPR shows how the size of the computational box and the presence 
of long waves affects the results on the pair-wise velocities. If we compare results for dark 
halos in the CHDM2 with those given by KHPR (bottom panel in Figure 15), then we find 
the difference of about 15% at l/i _1 Mpc. The smaller boxes used by KHPR give progressively 
lower estimates at larger radii, and the KHPR results for CDM were equally affected by the 
box size. If we take the same mass thereshold for dark halos, m\ > 1.3 X 10 10 h~ x Mq , the 
difference between CDM 6 = 1.5 simulations is about 20% at l/i _1 Mpc. We find the same 
difference for dark halos (about 20%). 

8. DISTRIBUTION OF DENSITY 

In CHDM, we of course have two different kinds of dark matter, cold and hot. Here 
we consider the distribution of each kind of particle as a function of total density. See 
also BHNPK, where the distributions are discussed further and visualized, both by showing 
isodensity surfaces for cold and hot densities and by showing the relative densities on cutting 
planes, for the CHDM2 simulation. 

The top panel in Figure 16 presents the ratio of the density of hot particles, /?hot, to the 
total density, /?total = Phot + Pcold, as a function of the total local density. The solid curve is 
for unsmoothed density and the dashed curve is for the density effectively smoothed with a 
Gaussian filter with a = 0.4 in units of our PM cell size (see below). It is interesting to note 
that cold and hot particles are distributed differently. In high density areas (/?total > (1 — 
2)/?total(/5total)) the density of hot particles is systematically lower by about 20% (phot / Ptotal ~ 
0.25) than the mean value = 0.30. Perhaps more surprising is the fact that the fraction 
of mass in hot particles does not change with density, for overdensity greater than about 10. 
(NKP95 shows that a related quantity, p co ld/ Phot, does rise near the centers of groups.) One 
might naively expect that the higher the density of an object, the larger is the fraction of 
neutrinos it captures because it has a deeper potential well. This is not what we see in our 
simulations. It seems that another effect plays a significant role: the larger is the density, the 
earlier it was formed and the larger was the velocity dispersion of neutrinos it accreted. At 
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the same time low density areas (lower than the mean density) have an excess of neutrinos. 
The enhancement is very large - a factor of two or more. It is readily apparent in the BHNPK 
visualizations. The much higher abundance of hot particles in low density regions no doubt 
impedes the gravitational collapse of galaxies and the development of larger scale structures 
in such regions. 

An immediate implication of this difference between hot and cold particles is that one 
must include hot and cold particles in numerical simulations in order to make accurate pre- 
dictions for CHDM models. It is not enough to have only the initial CHDM power spectrum, 
but evolve it with cold particles as in CDM models. Not only would the rate of evolution in 
the linear regime be wrong, but also the dynamics of both high and low density regions in the 
nonlinear regime would not be simulated properly. We note that in a number of publications 
(e.g., Mo et al. 1993; Croft k Efstathiou 1994; Gelb et al. 1993; and part of Jing et al. 1994) 
hot particles were not simulated and, thus, results of those publications do not accurately 
represent the CHDM model. 

The bottom panel of Figure 16 presents the density distribution function, P(p) } multi- 
plied by the total density p for the CHDM2 simulation (solid curve). For comparison we also 
present results from KHPR (dot-dashed curve). The agreement between the two simulations 
is very good in the density range p/ < p >= 10-1000. There is a lack of high density areas 
in KHPR simulations for pj < p » 1000, but the difference is not that large and is quite 
understandable taking into account the much better resolution in our present simulations. 
Effects of particle discreteness are seen as wiggles in the KHPR curve. 

9. CONCLUSIONS AND DISCUSSION 

In the present paper we have explained how we simulated a suite of cosmological models 
with = 1 and h = 0.5: Cold Dark Matter (CDM) with linear bias factor b = a^ 1 = 1 and 
1.5, and two simulations of Cold + Hot Dark Matter (CHDM) with tt c /tt„/tt b = 0.6/0.3/0.1, 
one started with the same random numbers as the CDM simulations (with large waves some- 
what higher than typical), the second with different random numbers corresponding to lower 
(more typical) amplitude for the largest waves. These simulations are a more detailed exam- 
ination of the models simulated in KHPR in a larger range of boxes but at lower resolution. 
Here we simulated periodic boxes of linear size 100 Mpc. We have analyzed these higher reso- 
lution (512 3 Particle Mesh) simulations both in the present paper and in a number of related 
papers, in particular NKP94, NKP95, and our visualization paper BHNPK; also Yepes et al. 

1994) , Ghigna et al. (1994), Bonometto et al. (1995), Dave et al. (1995), and Hellinger et al. 

1995) . 

In the present paper, we have mainly analyzed the evolution of the power spectrum of 
these models in real space, compared the power spectrum to observational data in redshift 
space, compared correlation functions of various populations in these simulations in both real 
and redshift space and with relevant data, analyzed the real-space rms velocity of pairs of 
galaxies cry in various simulations and compared to data; and analyzed the density of hot and 
cold particles as a function of total density. Our conclusions follows: 

Evolution of simulations. The evolution of our simulations is well described by linear 
theory on the largest scales, and by standard clustering analyses on the smallest scales 
we can resolve. 

Power spectra. We find that the power spectrum of the CHDM model in the nonlinear 
regime agrees rather well both with CfA2 data (if we compare this optical data with 
the mass-wighted power spectrum, to compensate for overmerging in dense regions) and 
with IRAS data (if we compare with the number-weighted power spectrum, since IRAS 
galaxies avoid dense regions). Moreover, with appropriate rescaling to account for the 
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biasing of these two different galaxy data sets, we find resonable agreement also in the 
linear regime. 

Correlation functions. Real-space autocorrelation functions £(r) were presented for 
CHDM dark matter halos ("galaxies") having masses greater than various lower limits. 
As expected, the higher-mass "galaxies" are more correlated. Linear and quasi-linear 
analyses appear to account for the behavior of £(r) on large scales. For CDM, £(r) for 
the dark matter is not well fit by a power law, but the dark halos were better fit by a 
power law and were anti-biased with respect to the dark matter. Redshift-space CHDM 
corrected on large scales by the method of Kaiser (1987), agree well with the CfA2 
data on all scales. But the CfA2 data does not agree as well with from the CDM 
simulation, £cdm(s) is too low especially on small scales s < 5/i _1 Mpc, and also on large 
scales s > 10/i _1 Mpc. 

Velocities. Here we concentrated on two statistics in real space: the rms velocity along 
a line joining a pair of halos <7||(i?), and rms velocities a v of dark matter and dark halos. 

We found &\\(R) to be systematically higher in CHDMi than CHDM2, reflecting the 
greater large-scale power in CHDMi, but lower in both CHDM simulations than in the 
unbiased CDM one. We see a fairly strong tendency for more massive pairs of "galaxies' 
to have larger pairwise velocities, although the rms velocities of the more massive halos 
is only slightly higher than that of the lower-mass halos. We measure a global velocity 
bias of about 0.8 for CDM and 1.0 for both CHDM simulations, but the pairwise velocity 
bias is 6 U)Pa irs = 0.44 — 0.62 in both CDM and CHDM simulations, with 6 U)Pa irs smaller 
for smaller-mass pairs of galaxies. 

Hot vs. Cold Particle Distribution. We find that the ratio of Hot/Total density is 
remarkably constant for overdensities exceeding about 10, and reduced compared to 
the average density in the box by about 20%. While hot particles are thus relatively 
depeleted in the dense regions, we find that their abundance is enhanced by a factor of 
two or more in the low-density regions. Thus the voids are hot dark matter balloons! 
It follows from these strikingly different distributions of cold and hot dark matter that 
accurately simulating CHDM requires treating the two different types of dark matter 
particles differently, as we have done. 

On the basis of the results discussed above, CHDM with fi c /fij, /fi& = 0.6/0.3/0.1 (i.e., 
with the hot dark matter consisting of neutrinos of mass 7 eV) appears to be a good fit to 
the data on the distribution and velocities of galaxies at the present epoch. We explained 
in NKP95, Ghigna et al. (1994), and especially in KBHP and Primack et al. (1995) why we 
think that a slightly lower neutrino mass ~ 5 eV may be preferred; but it remains to be seen 
whether the velocities on small scales will be too high, the void probability function will agree 
better with observations, etc. We are presently running and analyzing simulations with lower 
total neutrino mass. 
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FIGURE CAPTIONS 



Figure 1 Example of a large group of galaxies identified in the CHDM2 simulation. Coordi- 
nates are given in units of Mpc for h = 0.5. Positions of all cold particles inside a 3/i -1 Mpc 
cube are shown together with positions of dark halos with central overdensity larger than 
Sp/p > 150. Each circle on the plots corresponds to a dark halo, with the area of the circle 
being proportional to the mass of the dark halo. The group has mass M = 7.5 X 10 13 h~ x Mq , 
with nine "galaxies" exceeding the density threshold. The rms 3D velocity of cold particles 
in the group is 644 km s 1 . 

Figure 2 The same as Figure 1 , but for a small group of galaxies: the cold particle rms 3D 
velocity is 310 km s , M = 1.85 x 10 13 // '.\/. . The group has four members of about 
equal mass. Five "galaxies" on the upper part of the plot and one small dark halo at the 
bottom are not members of the group. 

Figure 3 The initial power spectra at z = 15. Because of the small number of harmonics 
with long waves, statistical fluctuations are apparent relative to theoretical spectra at k < 
lh Mpc" 1 . 

Figure 4 The evolution of the power spectrum in the CHDMi model. The solid curves 
show the power spectra of the total density ("cold" plus "hot") in the simulation at different 
redshifts. The dot-dashed curves present the linear theory predictions. Due to a statistical 
fluke, the first harmonic in the simulation (k = 0.15 h Mpc -1 ) is above the linear prediction 
just at the beginning of simulations. The dashed lines correspond to the power law slopes 
n = —1.2 and n = —2.5. 

Figure 5 Comparison of power spectra for CDM and CHDMi simulations in the nonlinear 
regime at the final moment (z = 0). The spectrum of the CHDM and that of the biased 
CDM (erg = 0.66) models are quite close in the range k = (0.1 — 1) h Mpc -1 , but the slope 

of the power spectrum in the nonlinear regime k = (0.2 — 2.0) h Mpc -1 is shallower for the 
CDM models (n —1) as compared with the CHDM spectrum, which results in a steeper 
correlation function at small radii for the CDM model. The dashed curve shows the linear 
spectrum for the biased CDM model. 

Figure 6 Comparison of the power spectrum of the galaxies in the CfA2 catalog (dots with 
error bars: Vogeley et al. 1992) with the spectrum predicted by the CHDM model in redshift 
space. The solid and dotted curves represent results for CHDMi and CHDM2 simulations 
respectively. The dashed curve is the linear power spectrum scaled upwards by factor 2.75 
(see text). 

Figure 7 The same as in Figure 6 , but for IRAS galaxies (dots are data of Fisher et al. 
1992, broken dashed curve are the same data as analyzed by Feldman et al. 1993). The linear 
CHDM spectrum (dot-dashed curve) was shifted upward by factor 1.47. 

Figure 8 The correlation functions £ of different populations in the CHDMi model at the 
present epoch (z = 0). The lower solid curve presents the correlation function of the "cold" 
dark matter in the simulation. In the range R = (0.5 — 5)/i -1 Mpc the correlation function 
follows the power law ^(R) = (i?/4.0/i -1 Mpc) -1 ' 8 indicated by the lower dot-dashed line. 
At larger radii the correlation function goes above the power law, but it remains reasonably 
close to the correlation function predicted by the linear theory for an infinite box (the dashed 
curve). The three upper solid curves on the plot show the correlation functions of dark 
halos chosen for different mass limits. From the bottom to the top, the curves are for the 
following mass limits and numbers of "galaxies": m\ > 1.28 X 10 10 h~ x Mq (corresponding 
to p > 50/o c ), N = 14606, m 1 > 2.56 x 10 10 /i -1 M , N = 6172, m 1 > 7.67 x 10 10 /i -1 M , 
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N = 1410. Here m\ indicates that the mass corresponds to the mass in the single cell at the 
density maximum. 

Figure 9 The correlation function of the "cold" dark matter in the CHDMi simulation is 
shown as the lower solid curve. The dashed curve represents the linear theory prediction. 
The same procedure was applied to the correlation function of dark halos with m\ > 1.28 X 
10 10 h~ x Mq . In this case the power spectrum predicted by the linear theory was multiplied 
by a factor 2.3 to account for the biasing of dark halos observed in Figure 8 . The dot-dashed 
lines on the plot are the same power-laws as in Figure 8 . 

Figure 10 The correlation function in real space for the unbiased CDM simulation. The 
correlation function of the dark matter deviates significantly from a power law, but that 
of the dark halos — which lies below that of the dark matter, i.e. it is "anti-biased" — 
does show power-law behavior (the dot-dashed line) with —1.8 slope and correlation length 
r ga _i = 5.7/i _1 Mpc. The same mass threshold (mi > 1.28 x 10 10 /i _1 M ; N = 22275) and 
weighting scheme was used as for the CHDM halos. No corrections were applied. 

Figure 11 The correlation function in the redshift space for the dark matter and dark halos 
in the CHDMi simulation. Note that at small radii the difference between halos and dark 
matter becomes larger. 

Figure 12 Comparison of predicted correlation function of "galaxies" in the CHDM model 
with the CfA2 data. The correlation functions in redshift space for the CHDMi model are 
shown as the solid curves. Corrections for the finite box size and the statistical fluctuations 
of few longest waves in the simulations were applied. The dashed curve represents the linear 
theory prediction amplified by the velocity effects. Triangles show the estimates of the corre- 
lation function for the CfA2 Redshift Survey (Vogeley et al. 1992). The agreement between 
the observational data and the model is reasonably good on all scales. 

Figure 13 The correlation function in redshift space for the unbiased CDM model. No 
corrections were applied. 

Figure 14 The correlation function in redshift space for the unbiased CDM model corrected 
for the box size and Kaiser's correction for 6 = 1. The dashed curve, representing linear 
theory in redshift space, does not match from the simulation very well. We tried many 
parameters, and this is the best match we could get. The triangles show the CfA2 data 
of Vogeley et al. (1992). The correlation function of the dark matter and the dark halos 
("galaxies") match very well at scales smaller than s < 3/i _1 Mpc, but £ ga i is about 1.5 times 
higher than ^ m for s < 3/i _1 Mpc. 

Figure 15 Real space rms velocity along the line joining a pair of objects. The top panel 
shows C||(-R) for dark halos with central overdensity larger than 100. The dot-dashed curves 
are for the CDM simulations (higher curve: 6 = 1.0, lower: 6 = 1.5). The bottom panel shows 
results for the CHDM2 simulation (solid curves). Different curves correspond to different 
mass thresholds: the larger the mass of both "galaxies" in the pair, the larger is their relative 
velocity. For comparison, we show the pair-wise rms velocity from the KHPR simulations for 
the CHDM model (dashed curve). 

Figure 16 Distribution of hot and cold particles in the CHDM2 simulation. The top panel 
shows the ratio of the density of hot particles to the total density as a function of the total 
local density. The bottom panel shows the density distribution function multiplied by the 
total density p for the CHDM2 simulation (solid curve) compared to the same quantity from 
KHPR. In order to make more reliable estimates for the low-density tail of P(p), we smoothed 
the total density field with a filter which weights the density of the central cell with weight 
3 and the six nearest cells with unit weight. The filter reduces discreteness effects due to 
small number of points in low density areas ("voids"). The dashed curve in the top panel 
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and both curves in the bottom panel in Figure 16 show results for the smoothed density; the 
unsmoothed density (obtained with the usual Cloud-In- Cell method) in shown as the solid 
curve in the top panel. 
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TABLE 1 

Velocity dispersions and Velocity biases 
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